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Abstract. Building on insights about the classical pathways to double ionization in 
strong laser fields we here propose a 1+1-dimensional model that captures essentials 
of the 3-d potential landscape and allows efficient studies of the process in reduced 
dimensions. The reduction to one degree of freedom for each electron is obtained by 
confining their motion to lines at an angle of 7r/6 with respect to the field axis; the 
justification for this choice is that upon variation of the electric field the Stark saddles 
move along the same lines. In this way we obtain a low-dimensional representation in 
which symmetric electron escape is possible. As a demonstration of the advantages of 
this model we confirm numerically the equivalent of the Wannier threshold behaviour 
for constant electric field and identify several classes of trajectories that can contribute 
to the ionization process. 

PACS numbers: 32.80.Rm, 32.80.Fb 

Introduction. Analyzing the dynamics of two-electron systems remains a 
challenge because of the high dimensionality of the configuration space, especially when 
highly excited states are involved, as in the case of strong field double ionization. Taylor 
group managed to obtain results for a full 3 +3- dimensional representation, albeit mostly 
at short wave length PUEJ- 1+1-dimensional models with both electrons aligned along 
one axis jSl El El E] are attractive (as reflected in a large number of publications based 
on this model), but physically flawed since the electron repulsion prevents an escape in 
the experimentally observed subspace of equal momenta [7J |HJ US EH HH EH EU • Among 
the alternatives that have been considered are models where the center of mass of the 
electrons is confined to move along the axis ^3]: they might be classified as 1.5+1.5 
dimensional, since the position of an electron along the field axis qualifies as a full 
degree of freedom and the perpendicular one, being shared with the other electron, as 
half a degree of freedom. It is our aim here to present and analyze a model with further 
constraints on the dynamics; the model nevertheless captures several essential features 
of the full dynamics and allows efficient simulations. 

The model is suggested by our previous analysis of the classical pathways to non- 
sequential double ionization in strong laser fields [131111 El- We found that the observed 
symmetric escape can be explained by the escape over a symmetrically arranged saddle. 
As the field changes, the saddle moves along lines that keep a constant angle with 
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respect to the polarization axis. Confining the electrons to run along 1-d tracks that 
pass through the saddles gives a model for 1+1 dimensional electron dynamics that has 
a potential landscape topology very similar to that for electrons in the full 3-d case. 
While the potential turns out to be very similar to the so called aligned-electron model 
[SI HI El E], it has one significant difference: since the electron tracks separate as they 
move away from the nucleus electron repulsion diminishes as the electrons go out. Thus, 
the diagonal where the two electron coordinates are equal, is accessible and not, as in the 
aligned model, suppressed by electron repulsion. This allows to mimic a key feature of 
the observed double ionization dynamics, namely the preference for symmetric electron 
escape [HI El EEI1 El EE21 EE! ■ 

In the present publication we analyze the model classically for the case of a static 
electric field. The static field forms a local maximum in the potential energy (that 
corresponds to the saddle in the 3-d case [THj) close to which simultaneous escape of 
electrons takes place. If the energy of the system equals the energy of the stationary 
point only purely symmetric electron motion leads to the simultaneous escape. For 
higher energy some deviations from the symmetric motion are allowed. The stability 
analysis of the saddle point in the full 3-d case allowed us to predict the dependence of the 
classical cross section for the process on energy close to the threshold jTTj, i.e. to obtain 
the Wannier threshold law QUI E0| m the presence of the external field. However, 
the results of the local analysis were difficult to verify in 3-d numerical simulations. The 
key features of the 3-d potential are present in our 1+1 dimensional model that allows 
us to obtain also the threshold law and moreover verify it in numerical simulations. 

The model. We begin with a reminder of the saddle configurations in the 
symmetric subspace in the 3-d case [IHl EHl E| ■ Assume the external electric field points 
in the ^-direction and the electrons are labeled i = 1,2 with positions r$ = (x^y^Zi). 
Then the Hamiltonian, in atomic units, is 

Let the electrons be moving in the x-z-plane, symmetric with respect to the field axis. 
Then their positions are (x, 0, z) and (— x, 0, z) and the potential energy becomes 

V = ; 4 + — j — r — 2Fz. (2) 

Vx 2 + z 2 2\x\ K ' 

This potential energy has a saddle at x s = 3 1 / 4 /(2\/F), z s = 3 3 / 4 /(2 v / F), with energy 
—3 3 / 4 2y/~F. Note that if we allow for variation of the field strength (like, e.g., in a laser 
field) the saddle moves along lines with z s /x s = \^3 =const. Because simultaneous 
electron escape takes place in the vicinity of the saddle, the idea is then to assume 
that the electrons move in a plane, yi = 0, with coordinates constrained exactly by this 
relation, zf = 3xj. This leaves only one degree of freedom for each electron, called r\ 
and r 2 : 

1^3. 1^3 
xi = -n z i = ~Y Tl and %2 = ~2 T2 Z2 = IT 7 * 2 ' 
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The Hamiltonian becomes 



pj + pj 2 2 1 Fy/3 

H = 7, |— | ~ |— | + / 7T~ ( r l+ r 2j- (4) 

* | r l| \ r 2\ J (n — r 2 ) 2 + T\T 2 



The Hamiltonian Q defines our 1+1 dimensional model and in following we restrict 
ourselves to description of the two electron system within this model. 

When the external field is present the potential energy in has a saddle located 

at 

3 i/4 

n = r 2 = r s = —= , (5) 



with energy 

V s = -3 3 / 4 2Vf. (6) 

The equipotential contours for F = 0.02 are plotted in Fig. ^ the stationary point near 
n = r2 ~ 9.3 is clearly visible. There is a symmetric subspace in the full phase space of 
the 1+1 dimensional system and the stationary point lives in this subspace. That is, if 
initial conditions are chosen symmetrically, i.e. t\ = r 2 and p\ = p 2 , the electron motion 
remains symmetric in time evolution. Trajectories that pass close to the stationary point 
and sufficiently close to the symmetric subspace lead to simultaneous electron escapes. 

Wannier threshold law. The potential energy in (j3J) possesses a stationary point 
located at (0), i.e. in the symmetric subspace of the system. The stability analysis of 
the stationary point in the full phase space allows us to determine the behaviour of the 
cross-section for simultaneous electron escape close to the threshold [T7| . There are two 
Lyapunov exponents of the point. The first, 

Ar = 3 l/8 j p3/4 ; (?) 
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with the eigenvector in the symmetric subspace, corresponds to a simultaneous motion 
of the electrons in the same direction away from the nucleus. Borrowing terminology 
from chemical reactions, we call this subspace the reaction coordinate. Because of the 
repulsion between electrons, there is an additional unstable direction, with Lyapunov 
exponent 

X± = ^3^F^ , (8) 

which enters in the threshold law. If the initial energy of the system equals the 
stationary point energy only a trajectory living in the symmetric subspace can lead 
to a simultaneous escape of both electrons. This reduces the dimensionality of the 
problem and the cross-section vanishes. For higher energy some deviations from the 
symmetric motion are possible, giving a finite volume of initial conditions and a non- 
vanishing cross-section. The dependence of the cross-section on energy o~(E) close to 
the stationary point energy V s can be obtained in the spirit of the Wannier analysis 
Pll [El EH EDI EQ , resulting in 

a{E) ~(E- V s y, (9) 

with an exponent 

"=aHtT~ L354 ' (10) 

This is larger than the corresponding exponent in the full 3-d case where, for a doubly 
charged remaining ion, the exponent is 1.292 The cross section is larger if the 

exponent is smaller, i.e. if the saddle is crossed more quickly (larger A r ) or if the 
differences from the symmetric motion grow more slowly (smaller Aj_). This cross section 
exponent is an additional characteristic of the double ionization process. 

Numerical results. The derivation of the cross section exponent (fT0|) for the 
simultaneous electron escape is based on the local analysis of the stationary point of the 
system. One may wonder if the higher order terms can modify the behaviour of the cross 
section. To test the results of the local analysis we can perform numerical simulations 
of the process. In the 3-d case considered in Ref. [T7] that was quite difficult due to 
high dimensionality of the problem. Here, we deal effectively with the three dimensional 
phase space and the numerical simulations become feasible. 

In order to avoid problems with the Coulomb singularities, we can add cut-offs in 
the denominator, so that the potential becomes 

V = -- j=L= - -=L= + 1 - ^(n + r 2 ).(ll) 

yVi + 1 yrl + l \J (n - r 2 ) 2 + nr 2 + 1 1 

The cut-offs change slightly the prediction for the threshold exponent which now equals 
fi = 1.357 (for F = 0.02). We have run trajectories with initial conditions chosen 
microcanonically for different energies E above the threshold energy V s but with an 
additional requirement that they have to lie on the surface r\ + r 2 = 0. We deal with 
an open system and the condition t\ + r 2 = ensures that the trajectories start in 
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Figure 2. Number of ionized trajectories (circles with statistical error bars) versus 
energy above the threshold energy, i.e. the energy of the maximum of the potential l|llfl , 
corresponding to F = 0.02 (1.5 x 10 6 trajectories have been integrated to obtain each 
point). Solid line is the fit of © to the numerical data that results in fi = 1.383±0.031. 

the vicinity of the nucleus. In Fig. El number of trajectories leading to simultaneous 
electron escape together with the numerical fit of the function Q versus energy are 
shown. The fitted value of the cross section exponent is 1.383 ± 0.031 what agrees with 
the theoretical prediction. 

In Fig. El we have included only trajectories leading to simultaneous escape of both 
electron. However, it is possible that the first electron ionizes by passing close to the 
stationary point while the other one is returned and after revisiting the nucleus ionizes 
too. Actually there is a separatrix that divides the phase space into areas corresponding 
to simultaneous escapes and escapes with a single revisit of the nucleus. The separatrix 
consists of trajectories where the first electron escapes while the other one approaches a 
single electron Stark saddle and looses its kinetic energy. Finally, when the first electron 
is gone, the other electron remains at the single electron saddle with no kinetic energy. 
In Fig. El there are examples of ionized trajectories and trajectories belonging to the 
separatrix. All trajectories in Fig. El correspond to: F = 0.02, energy E = V s + 0.1, 
symmetric initial momenta pi = P2 and positions lying on a circle, i.e. r\ = {r s /2) cos a 
and r\ = {r s /2) sina. With these conditions if | tt/4 — a\ is smaller than about 0.01547T 
simultaneous escape takes place. For greater 1 7r/4 — a\ but smaller than about 0.01557T 
the second electron revisits the nucleus and then escapes. At |7r/4 — a\ ~ 0.01557T 
another separatrix appears that separates double ionization trajectories with single and 
double revisit of the nucleus. This separatrix (see Fig. EJ) consists of trajectories where 
the second electron after revisiting the nucleus approaches the single electron Stark 
saddle loosing its kinetic energy so that finally it stands on the saddle. With further 
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Figure 3. Panel (a): trajectories (dashed lines) belonging to the separatrix that 
divides the phase space into areas of simultaneous electron escape and double escape 
with a single revisit of the nucleus by one of the electrons; solid lines are examples of 
trajectories leading to simultaneous electron escape. Panel (b): trajectories (dashed 
lines) belonging to the separatrix that divides double ionization trajectories with single 
and double revisit of the nucleus; solid lines are examples of trajectories leading to 
double ionization with a single revisit of the nucleus by one of the electrons. All 
trajectories correspond to: F = 0.02, energy E = V s + 0.1, symmetric initial momenta 
Pi = p2 and positions lying on a circle, i.e. r\ = (r s /2) cosa and r% = (r s /2) sina. 



increase of 1 7r/4 — a\, number of revisits increases quickly until the returning electron 
has too small energy to cross the single electron saddle even if the first electron is gone. 

Conclusions. The analysis of pathways to non-sequential double ionization of 
atoms in strong fields (THIIIEIIIZI allows us to propose a 1+1 dimensional model of the 
process. In the present publication we analyze the model in the case of a static electric 
field. The Wannier threshold law for simultaneous electron escape in the presence of 
the static field is derived and tested numerically. We also identify separatrices that 
divide the phase space of the system into areas of simultaneous electron escapes and 
double escapes with a multiple revisit of the nucleus by one of the electrons. In classical 
mechanics these different trajectories contribute independently to the double ionization 
process. In quantum mechanics one may expect interesting coherence effects resulting 
from interference of the different paths. 

The model is considered classically and for a static external field only. However, 
quantum calculations of double ionization in strong laser field |22[ EHl 121] including: 
electron tunneling, rescattering and subsequent double escape can be performed 
numerically very efficiently within the model. As we have mentioned in the introduction 
the crucial advantage of the model over the well known aligned electron model [31I11EIE] 
is that it does not forbid symmetric simultaneous escapes of electrons what is observed 
experimentally and what can not be described by the aligned electron model due to an 
overestimation of the Coulomb repulsion. 
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